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Abstract. We report the detection of a significant (5.5 a) 
excess of correlations between galaxy ellipticities at scales 
ranging from 0.5 to 3.5 arc-minutes. This detection of a 
gravitational lensing signal by large-scale structure was 
made using a composite high quality imaging survey of 
6300 arcmin 2 obtained at the Canada France Hawaii Tele- 
scope (CFHT) with the UH8K and CFH12K panoramic 
CCD cameras. The amplitude of the excess correlation is 
2.2 ± 0.2 % at 1 arcmin scale, in agreement with theoreti- 
cal predictions of the lensing effect induced by large-scale 
structure. We provide a quantitative analysis of systemat- 
ics which could contribute to the signal and show that the 
net effect is small and can be corrected for. In particular, 
we show that the spurious excess of correlations caused 
by the residual of the anisotropic Point Spread Function 
(PSF) correction is well below the measured signal. We 
show that the measured ellipticity correlations behave as 
expected for a gravitational shear signal. The relatively 
small size of our survey precludes tight constraints on cos- 
mological models. However the data are in favor of clus- 
ter normalized cosmological models, and marginally reject 
Cold Dark Matter models with (ft = 0.3, ct 8 < 0.6) or 
(fi = 1, (78 = 1)- The detection of cosmic shear demon- 
strates the technical feasibility of using weak lensing sur- 
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veys to measure dark matter clustering and the potential 
for cosmological parameter measurements, in particular 
with upcoming wide field CCD cameras. 

Key words: Cosmology: theory, dark matter, gravita- 
tional lenses, large-scale structure of the universe 



1. Introduction 

The measurement of weak gravitational lensing produced 
by the large-scale structures in the universe (hereafter, the 
cosmic shear) is potentially the most effective, albeit chal- 
lenging, step toward a direct mapping of the dark matter 
distribution in the universe at intermediate and low red- 
shift. Unlike several popular probes of large-scale struc- 
tures, lensing maps the dark matter directly, regardless of 
the distribution of light emitted by gas and galaxies or the 
dynamical stage of the structures analysed. 

A decade of theoretical and technical studies has 
shown that the gravitational distortion produced by the 
structures along the lines-of-sight contains important 
clues on structure formation models (see Mellier 1999, 
Bartelmann & Schneider 2000 for reviews and references 
therein). From these studies, we know that weak lens- 
ing can provide measurements of cosmological param- 
eters and the shape of the projected density power 
spectrum (Blandford et al. 1991, Miralda-Escude 1991, 
Kaiser 1992, Villumsen 1996, Bernardeau et al. 1997, 
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Jain & Seljak 1997, Kaiser 1998, Schneider et al. 1998, 
Jain et al. 2000, Van Waerbeke et al. 1999, 

Bartelmann & Schneider 1999). However, it is also clear 
that the most challenging issues are observationals, be- 
cause the measurement of extremely weak gravitational 
distortions is severely affected by various sources of noise 
and systematics such as the photon noise, the opti- 
cal distortion of astronomical telescopes and the atmo- 
spheric distortion. Therefore, the problem of reliable shape 
measurement has also received much attention in the 
last few years (Bonnet & Mellier 1995, Kaiser et al. 1995, 
Van Waerbeke et al. 1997, 

Hoekstra et al. 1998, Kuijkcn 1999, Rhodes et al. 1999, 
Kaiser 1999, Bertin 2000). 

Despite considerable difficulties in recovering weak 
lcnsing signals, the potential cosmological impact of the 
cosmic shear analysis has motivated several teams to de- 
vote efforts on imaging surveys designed for the mea- 
surement of the galaxy distortion produced by grav- 
itational lensing, either by observing many indepen- 
dent small fields, like the VLT/FORS-I (Maoli et al in 
preparation), the HST/STIS (Seitz et al 1998), the WHT 
(Bacon et al. 2000), or by observing few intermediate to 
large fields, like the CFHT / CFH12K-UH8K (this work 
and Kaiser et al. 2000), the SDSS (Annis et al 1998) and 
other ongoing surveys. In this paper we present the re- 
sults of the analysis based on 2 square degrees obtained 
during previous independent observing runs at the CFHT 
with mixed I and V colors. This study is part of a our 
weak lensing survey carried out at CFHT (hereafter the 
DESCART project 1 ) which will cover 16 square degrees 
in four colors with the CFH12K camera. Though the sur- 
vey is far from completion, data obtained during previous 
runs have been used jointly with the first observations of 
the DESCART survey that we did in May 1999 and in 
November 1999 in order to demonstrate that the techni- 
cal issues can be overcome and to better prepare the next 
observations. This set of data permits us already to report 
on the detection of a cosmic shear signal. 
In the following, we discuss the technique used to extract 
the cosmological signal and to measure its amplitude and 
show that systematic effects are well under control. The 
paper is organized as follows. Section 2 describes our data 
sets. Section 3 discusses the details of our PSF correction 
procedure and Section 4 presents the final results. Section 
5 is devoted to the discussion of the residual systemat- 
ics and their correction. Section 6 presents a preliminary 
quantitative comparison of our signal with numerical ex- 
pectations of cosmological scenarios as derived from ray- 
tracing simulations. Conclusions are given in Section 7. 



2. Description of the data 

The difficulty to get a wide angle coverage of the sky in 
good conditions is the reason why there is not yet a clear 
detection of cosmic shear. For this work, we decided to get 
the widest angular field possible, which was done at the 
expense of homogeneity of the data set. However this does 
not impact our primary goals which are the detection of a 
weak lensing signal and the test of the control of system- 
atics. 

We use in total eight different pointings mixing CFH12K 
and UH8K data sets (see Table 1). They are spread 
over five statistically independent areas, each separated 
by more than 10 degrees. The total field covers about 
6300arcmin 2 , and contains 3 x 10 5 galaxies (with a number 
density n g ~ 30 gal/arcmin 2 ). Note that the galaxies are 
weighted as discussed in Section 2.2, and parts of the fields 
are masked, so the effective number density of galaxies is 
about half. 

All the data were obtained at the CFHT prime fo- 
cus. We used observations spread over 4 years from 
1996 to 1999, with two different cameras: the UH8K 
(Luppino et al 1994), covering a field of 28x28 square arc 
minutes with 0.2 arc-second per pixel and the CFHT12K 
2 (Cuillandre et al 2000) covering a field of 42x28 square 
arc-minutes with 0.2 arc-second per pixel as well. Because 
these observations were initially done for various scientific 
purposes, they have been done cither in I or in V band. Ta- 
ble 1 summarizes the dataset. The SA57 field was kindly 
provided by M. Creze and A. Robin who observed this 
field for another scientific purpose (star counts and proper 
motions). The UH8K Abell 1942 data were obtained dur- 
ing discretionary time. The F14 and F02 fields are part of 
the deep imaging survey of 16 square-degrees in BVRI be- 
ing conducted at CFHT jointly by several French teams. 
This survey is designed to satisfy several scientific pro- 
grams, including the DESCART weak lensing program, 
the study of galaxy evolution and clustering evolution, 
clusters and AGN searches, and prepare the spectroscopic 
sample to be studied for the VLT-VIRMOS deep red- 
shift survey (Le Fevre et al 1998). CFDF-03 is one of the 
Canada-France-Deep-Fields (CFDF) studied within the 
framework of the Canada-France Deep Fields, with data 
collected with the UH8K (Mc Cracken et al in prepara- 
tion). 

The observations were done as usual, by splitting the to- 
tal integration time in individual exposures of 10 minutes 
each, offsetting the telescope by 7 to 12 arc-seconds after 
each image acquisition. For the I and the V band data, we 
got between 7 to 13 different exposures per field, all with 
seeing conditions varying by less than ± 0.07 arc-seconds 
(the others were not co-added). The total exposure times 
range from 1.75 hours in V to 5 hours in I. 
The total field observed covers 2.05 square degrees, includ- 
ing 0.88 square degrees in V and 1.17 square degrees in I. 



1 http://terapix.iap.fr/Descart/ 
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Table 1. List of the fields. Most of the exposures were taken in the I band at CFHT. The total area is 1.7 deg 2 , and the 8 fields 
are uncorrclated. 



Target 


Name 


Camera 


Used area 


Filter 


Exp. time 


Period 


seeing 


F14P1 


Fl 


CFH12K 


764 arcmin^ 


V 


5400 sec. 


May 1999 


0.9" 


F14P2 


F2 


CFH12K 


764 arcmin 2 


V 


5400 sec. 


May 1999 


0.9" 


F14P3 


F3 


CFH12K 


764 arcmin 2 


V 


5400 sec. 


May 1999 


0.9" 


CFDF-03 


F4 


UH8K 


669 arcmin 2 


I 


17000 sec. 


Dec. 1996 


0.75" 


SA57 


F5 


UH8K 


669 arcmin 2 


I 


12000 sec. 


May 1998 


0.75" 


A 1942 


F6 


UH8K 


573 arcmin 2 


I 


10800 sec. 


May 1998 


0.75" 


F02P1 


F7 


CFH12K 


1050 arcmin 2 


I 


9360 sec. 


Nov. 1999 


0.8" 


F02P4 


F8 


CFH12K 


1050 arcmin 2 


I 


7200 sec. 


Nov. 1999 


0.9" 



However, one CCD of the UH8K and two CCDs initially 
mounted on the CFH12K of the May 1999 run have strong 
charge transfer efficiency problems and are not suitable 
for weak lensing analysis. Therefore, the final area only 
covers 1.74 square degrees: 0.64 square degrees in V and 
1.1 square degrees in I. As we can see from Table 1 each 
field has different properties (filter, exposure time, seeing) 
which makes this first data set somewhat heterogeneous. 
The data processing was done at the TERAPIX data cen- 
ter located at IAP which has been created in order to 
process big images obtained with these panoramic CCD 
cameras 3 . Its CPU (2 COMPAQ XP1000 with 1.2 Gb 
RAM memory each equipped with DEC alpha ev6/ev67 
processors) and disk space (1.2 Tbytes) facilities permit 
us to handle such a huge amount of data efficiently. 
For all but the CFDF-03 field, the preparation of the de- 
trending frames (master bias, master dark, master flats, 
superflats, fringing pattern, if any) and the generation of 
pre-reduced and stacked data were done using the FLIPS 
pre-reduction package (FITS Large Image Pre-reduction 
software) implemented at CFHT and in the TERAPIX 
pipeline (Cuillandre et al in preparation). In total, more 
than 300 Gbytes of data have been processed for this work. 
The CFHT prime focus wide-field corrector intro- 
duces a large-scale geometrical distortion in the field 
(Cuillandre et al 1996). Re-sampling the data over the an- 
gular size of one CCD (14 arc-minutes) cannot be avoided 
if large angular offsets (> 40 arc-seconds) are used for the 
dithering pattern (like for the CFDF-03 data). Since we 
kept the offsets between all individual exposures within a 
15 arc-seconds diameter disk, the contribution of the dis- 
tortion between objects at the top and at the bottom of 
the CCD between dithered exposures is kept below one 
tenth of a pixel. With the seeing above 0.7 arc-second and 
a sampling of 0.2 arc-second/pixel, the contribution of this 
effect is totally negligible. A simulation of the optical dis- 
tortion of the instrument shows that the variation from 
one field to another never exceeds 0.3%, which confirms 
what we expected from the CFHT optical design of the 
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wide-field corrector. We discuss this point in Section 5, in 
particular by confirming that the sensitivity of the shear 
components with radial distances is negligible. Also not 
correcting this optical distortion results in a slightly dif- 
ferent plate scale from the center to the edge of the field 
(pixels see more sky in the outside field). But this is also 
of no consequence for our program since the effect is very 
small as compared to the signal we are interested in. 
The stacking of the non-CFDF images has been done in- 
dependently for each individual CCD (each covering 7x14 
arc- minutes). We decided not to create a single large 
UH8K or CFH12K image per pointing since it is useless 
for our purpose. It complicates the weak lensing analysis, 
in particular for the PSF correction, and needs to handle 
properly the gaps between CCDs which potentially could 
produce discontinuities in the properties of the field. The 
drawback is that we restricted ourselves to a weak lens- 
ing analysis on scales smaller than 7 arc-minutes (radius 
smaller than 3.5 arc-minutes, as shown in the next fig- 
ures); this is not a critical scientific issue since the total 
field of view is still too small to provide significant signal 
beyond that angular scale. In the following we consider 
each individual CCD as one unit of the data set. 
The co-addition was performed by computing first the off- 
set of each CCD between each individual exposure from 
the identification of common bright objects (usually 20 ob- 
jects) spread over one of the CCD's arbitrary chosen as a 
reference frame. Then, for each exposure the offsets in the 
x- and y- directions are computed using the detection algo- 
rithm of the SExtractor package (Bertin & Arnouts 1996) 
which provides a typical accuracy better than one tenth 
of a pixel for bright objects. The internal accuracy of this 
technique is given by the rms fluctuations of the offsets of 
each reference object. Because our offsets were small the 
procedure works very well and provides a stable solution 
quickly. We usually reach an accuracy over the CCDs of 
0.25 pixels rms (0.05 arc-second) in both directions for 
offsets of about 10 arc-seconds (50 pixels). Once the off- 
sets are known the individual CCDs are stacked using a 
bilinear interpolation and by oversampling each pixel by 
a factor of 5 in both x- and y- directions (corresponding 
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to the rms accuracy of the offsets). The images are then 
re-binned lxl and finally a clipped median procedure is 
used for the addition. The procedure requires CPU and 
disk space but works very well, provided the shift between 
exposures remains small. We then end up with a final set 
of stacked CCDs which are ready for weak lensing analy- 
sis. 

The twelve separate pointings of the CFDF-03 field were 
processed independently using a method which is fully de- 
scribed elsewhere (McCracken et al 2000, in preparation). 
Briefly, it uses astrometric sources present in the field to 
derive a world coordinate system (WCS; in this work we 
use a gnomic projection with higher order terms). This 
mapping is then used to combine the eight CCD frames to 
produce a single image in which a uniform pixel scale is re- 
stored across the field. Subsequent pointings are registered 
to this initial WCS by using a large number of sources dis- 
tributed over the eight CCDs to correct for telescope flex- 
ure and atmospheric refraction. For each pointing the reg- 
istration accuracy is ~ 0.05" rms over the entire field. The 
final twelve projected images are combined using a clipped 
median, which, although sub-optimal in S/N terms, pro- 
vides the best rejection for cosmic rays and other transient 
events for small numbers of input images. 

3. Galaxy shape analysis 

The galaxies have been processed using the IMCAT soft- 
ware generously made available by Nick Kaiser 4 . Some of 
the process steps have been modified from the original IM- 
CAT version in order to comply with our specific needs. 
These modifications are described now. 

The object detection, cen- 

troid, size and magnitude measurements are done using 
SExtractor (Bcrtin & Arnouts 1996 5 ) which is optimized 
for the detection of galaxies. We replace the parameter 
r g (physical size of an object), calculated in the IMCAT 
peak finder algorithm by the half-light-radius of SExtrac- 
tor (which is very similar to r^ measured in IMCAT). 
This lowers the signal-to-noise of the shape measurements 
slightly, but it is not a serious issue for the statistical de- 
tection of cosmic shear described in this work. Before go- 
ing into the details of the shape analysis, we first briefly 
review how IMCAT measures shapes and corrects the stel- 
lar anisotropy. Technical details and proofs can be found in 
Kaiser ct al. 1995 (hereafter KSB), Hoekstra et al. 1998 
and Bartelmann & Schneider 2000. 

3.1. PSF correction: the principle 

KSB derived how a gravitational shear and an anisotropic 
PSF affect the shape of a galaxy. Their derivation is a 
first order effect calculation, which has the nice property 
to separate the gravitational shear and the atmospheric 
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effects. The correction is calculated first on the second 
moments of a galaxy, and subsequently the galaxy cllip- 
ticity can be directly expressed as a function of the shear 
and the star anisotropy. The raw ellipticity e of an object 
is the quantity measured from the second moments If,- of 
the surface brightness /(<?): 
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The aim of the window function W(0) is to suppress the 
photon noise which dominates the objects profile at large 
radii. According to KSB, in the presence of a shear 7^ and 
a PSF anisotropy pp, the raw ellipticity e is sheared and 
smeared, and modified by the quantity <5e: 



8e a = PZfo + PZ? Pfl . 



(2) 



The shear and smear polarization tensors P^p and P^™ 
are measured from the data, and the stellar ellipticity p, 
also measured from the data, is given by the raw stellar 
ellipticity e*: 



p a = 



psm 

aa 



(3) 



Using Eq.(2) and Eq.(3) we can therefore correct for the 
stellar anisotropy, and obtain an unbiased estimate of 
the orientation of the shear 7,3. To get the right ampli- 
tude of the shear, a piece is still missing: the isotropic 
correction, caused by the filter W{6) and the isotropic 
part of the PSF, which tend to circularize the objects. 
Luppino & Kaiser 1997 absorbed this isotropic correction 
by replacing the shear polarization P sh in Eq(2) (which is 
an exact derivation in the case of a Gaussian PSF) by the 
pre-seeing shear polarisability P 7 : 



jsh 



(4) 



This factor 'rescales' the galaxy ellipticity to its true 
value without changing its orientation, after the stellar 
anisotropy term was removed. The residual anisotropy left 
afterwards is the cosmic shear 73, therefore the observed 
ellipticity can be written as the sum of a 'source' elliptic- 
ity, a gravitational shear 7 term, and a stellar anisotropy 
contribution: 



gObs ^source 



+ P 7 7 + P sm p. 



(5) 



There is no reason that e source should be the 
true source ellipticity e true , as demonstrated by 
Bartelmann & Schneider 2000. The only thing we know 
about e source is that (e true ) = implies ( e source ) = 0. 
Therefore Eq.(5) provides an unbiased estimate of the 
shear 7 as long as the intrinsic ellipticities of the galaxies 
are uncorrelated (which leads to (e true ) = 0). The esti- 
mate of the shear is simply given by 
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= P" 1 • (e obs - P sm p). 



(6) 
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The quantities P 7 , P sm and p are calculated for each 
object. The shear estimate per galaxy (Eq(6)) is done us- 
ing the matrices of the different polarization tensors, and 
not their traces (which corresponds to a scalar correc- 
tion) as often done in the literature. Although the dif- 
ference between tensor and scalar correction is small (be- 
cause P 7 is nearly proportional to the identity matrix), 
we show elsewhere, in a comprehensive simulation pa- 
per (Erben et al. 2000), that the tensor correction gives 
slightly better results. 

3.2. PSF correction: the method 

The process of galaxy detection and shape correction can 
be done automatically, provided we first have a sam- 
ple of stars representative of the PSF. However, in prac- 
tice the star selection needs careful attention and cannot 
be automated because of contaminations. Stars can have 
very close neighbor(s) (for instance a small galaxy exactly 
aligned with it) that their shape parameters are strongly 
affected. Therefore we adopted a slow but well-controlled 
manual star selection process: on each CCD, the stars are 
first selected in the stellar branch of the r g — mag diagram 
in order to be certain to eliminate saturated and very faint 
stars. We then perform a 3a clipping on the corrected star 
ellipticitics, which removes most of the stars whose shape 
is affected by neighbors (they behave as outliers compared 
to the surounding stars). It is worth noting that the a clip- 
ping should be done on the corrected ellipticities and not 
on the raw ellipticities, since only the corrected ellipticities 
are supposed to have a vanishing anisotropy. The stellar 
outliers which survived the a clipping are checked by eye 
individually to make sure that no unusual systematics are 
present. 

During this procedure, we also manually mask the regions 
of the CCD which could potentially produce artificial sig- 
nal. This includes for example the areas with very strong 
gradient of the sky background, like around bright stars or 
bright /extended galaxies, but also spikes produced along 
the diffraction image of the spider supporting the sec- 
ondary mirror, columns containing light from saturated 
stars, CCD columns with bad charge transfer efficiency, 
residuals from transient events like asteroids which cross 
the CCD during the exposure and finally all the bound- 
aries of each CCD. At the end, we are left with a raw 
galaxy catalogue and a star catalogue free of spurious ob- 
jects, and each CCD chip has been checked individually. 
This masking process removes about 15% of the CCD area 
and the selection itself leaves about 30 to 100 usable stars 
per CCD. 

The most difficult step in the PSF correction is Eq.(6), 
where the inverse of a noisy matrix P 7 is involved. If we do 
not pay attention to this problem, we obtain corrected el- 
lipticities which can be very large and/or negative, which 
would force us to apply severe cuts on the final catalogue 
to remove aberrant corrections, thus losing many objects. 



A natural way to solve the problem is to smooth the ma- 
trix P 7 before it is inverted. In principle P 7 should be 
smoothed in the largest possible parameter space defining 
the objects: P 7 might depend on the magnitude, the ellip- 
ticity, the profile, the size, etc... In practice, it is common 
to smooth P 7 according to the magnitude and the size 
(see for instance Kaiser et al. 1998, Hoekstra et al. 1998). 
Smoothing performed on a regular grid is generally not 
optimal, and instead, we calculate the smoothed P 7 for 
each galaxy from its nearest neighbors in the objects pa- 
rameter space (this has the advantage of finding locally 
the optimal mesh size for grid smoothing) . Increasing the 
parameter space for smoothing does not lead to significant 
improvement in the correction (which is confirmed by our 
simulations in Erben et al. 2000), therefore we keep the 
magnitude and the size r^ to be the main functional de- 
pendencies of P 7 . 

A smoothed P 7 does not eliminate all abnormal ellip- 
ticities; the next step is to weight the galaxies according to 
the noise level of the ellipticity correction. Again, this can 
be done in the gridded magnitude/size parameter space 
where each cell contains a fixed number of objects (the 
nearest neighbors method). We then calculate the vari- 
ance a 2 of the ellipticity of those galaxies, which gives 
an indication of the dispersion of the ellipticities of the 
objects in the cell: the larger a~ , the larger the noise. We 
then calculate a weight w for each galaxy, which is directly 
given by a 2 : 

exp(— 5(cr £ — a) 2 ) if a e < 1 
W = I \ exp(-5(l - a) 2 ) if a e > 1 ' (7) 

where a is a free parameter, which is chosen to be the max- 
imum of the ellipticity distribution of the galaxies. Eq.(7) 
might seem arbitrary compared to the usual 1/of weight- 
ing, but the inverse square weighting tends to diverge for 
low-noise objects (because such objects have a small a 2 ), 
which create a strong unbalance among low noise objects. 
The aim of the exponential cut-off as defined in Eq.(7) is 
to suppress this divergence 6 . 

The weighting function prevents the use of an arbi- 
trary and sharp cut to remove the bad objects. However, 
we found in our simulations (Erben et al. 2000) that we 
should remove objects smaller than the seeing size, since 
they carry very little lensing information, and the PSF 
convolution is likely to dominate the shear amplitude. Our 
final catalogue contains about 191000 galaxies, of which 
23000 are masked. It is a galaxy number density of about 
n c± 26 gal/arcmin 2 , although the effective number density 
when the weighting is considered should be much less. We 
find a = 0.5, which corresponds to the ellipticity variance 
of the whole catalogue. 

6 Note that the use of a different weighting scheme like w oc 
1/(q 2 +(Je) has almost no effect on the detection. Other weight- 
ing schemes have been used, such as in Hoekstra et al. 2000 
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4. Measured signal 

The quantity directly accessible from the galaxy shapes 
and related to the cosmological model is the variance 
of the shear (j 2 ). An analytical estimate of it us- 
ing a simplified cosmological model (power-law power 
spectrum, sources at a single redshift plane, lead- 
ing order of the perturbation theory, and no cos- 
mological constant) gives (Kaiser 1992, Villumsen 1996, 
Bernardeau et al. 1997, Jain & Seljak 1997): 



From the unweighted galaxy ellipticities e a , an esti- 
mate of 7 2 (^i) at the position 8i is given by: 



a=l,2 \ k=l / 



(9) 



The inner summation is performed over the iV galaxies 
located inside the smoothing window centered on 9i, and 
the outer summation over the ellipticity components. The 
ensemble average of Eq.(9) is 



0.75 _0. 75 



larcmin 



<£[7 2 (0,)]H^ + <7 2 > 



(10) 



(8) 



where n is the slope of the power spectrum, as its nor- 
malization, z s the redshift of the sources and 9 the top- 
hat smoothing filter radius. The expected effect is at the 
percent level, but at small scales the non-linear dynam- 
ics is expected to increase the signal by a factor of a few 
(Jain & Seljak 1997). Nevertheless Eq.(8) has the advan- 
tage of clearly giving the cosmological dependence of the 
variance of the shear. 
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Fig. 1. Square-root of the variance of the measured shear 
as a function of the radius of the top-hat window (solid 
line). The maximum angular scale, 3.5 arc- minutes radius, 
is fixed by the maximum angular scale defined by individ- 
ual CCDs (7'). Error bars are computed over 1000 random 
realizations of the galaxy catalogue. The other lines are the- 
oretical predictions of the same quantity for different cosmo- 
logical models in the non-linear regime (using the fitting for- 
mula in Peacock & Dodds 1996): the long-dashed line corre- 
sponds to (Q = 1,A = 0,O8 = 0.6), the dashed line to 
(fl = 0.3, A = 0,(78 = 0.6), and the dot-dashed line to 
(fi = 0.3, A = 0.7,(78 = 0.6). 




The term <J 2 /N can be easily removed using a random 
realization of the galaxy catalogue: each position angle of 
the galaxies is randomized, and the variance of the shear 
is calculated again. This randomization allows us to de- 
termine <J 2 /N and the error bars associated with the noise 
due to the intrinsic ellipticity distribution. At least 1000 
random realizations are required in order to have a precise 
estimate of the error bars. Note that it is strictly equiv- 
alent to use an estimator where the diagonal terms are 
removed in the sum (9), which suppress automatically the 
a 2 /N bias. 

When we take into account the weighting scheme for 
each galaxy, the estimator Eq.(9) has to be modified ac- 
cordingly as follows: 



E\i 2 {e t )] = x: 

(1 = 1,2 



2_^w(9 k )e a (6 k 

fe=i 



N 



V fe=i 



(11) 



/ 



where w is the weight as defined in Eq.(7). The variance of 
the shear is not only the easiest quantity to measure, but it 
is also fairly weakly sensitive to the systematics provided 
that they are smaller than the signal. The reason is that 
any spurious alignment of the galaxies, in addition to the 
gravitational effect, adds quadratically to the signal and 
not linearly: 



ilmes) (itrue) ilbias) ' 



(12) 



Therefore, a systematic of say 1% for a signal of 3% only 
contributes to ~ 5% in (7 2 ) 1 / 2 . We investigate in detail in 
the next Sections the term ("f 2 ias ) and show that it has a 
negligible contribution. 

We will present results on the shear variance measured 
from the data sets described in Section 2. The variance 
ilmes) i s measured in apertures which are placed on a 
10 x 20 grid for each of the 2000 x 4000 CCDs. By construc- 
tion the apertures never cross the CCD boundaries, and 
if more than 10% of the included objects turns out to be 
masked objects, this aperture is not used. Figure 1 shows 
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Fig. 2. For different smoothing sizes (indicated at the top of each panel), the value of the measured signal (given by the arrow) 
compared to the signal measured in the randomized catalogues (histograms) . This figure shows how far the signal deviates from 
a pure random orientation of the galaxies. Note that the distribution of (j 2 ) is not Gaussian. 
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(imes) 1 ^ 2 (thick line) with error bars obtained from 1000 
random realizations. The three other thin lines correspond 
to theoretical predictions obtained from an exact numer- 
ical computation for three different cosmological models, 
in the non-linear regime. We assumed a normalized broad 
source redshift distribution given by 



n(z s ) 




wit t the parameters (zo,a,(3) — (0.9,2,1.5) are sup- 
posed to match roughly the redshift distribution in our 
data sets 7 . The shape of this redshift distribution mimics 
those observed in spectroscopic magnitude-limited sam- 
ples as well as those inferred from theoretical predictions 
of galaxy evolution models. Since we did significant selec- 
tions in our galaxy catalog the final redshift distribution 
could be modified. We have not quantified this, but we do 
not think it would significantly change the average red- 
shift of the sample, even if the shape may be modified. 
The variance of the shear (-j 2 ) is computed via the for- 
mula (see Schneider et al. 1998 for the notations): 

/>oo 

< 7 2 )= 27r/ k&kP K {k)I 2 TH {k8), (14) 
Jo 

where is the Fourier transform of a Top-Hat window 
function, and P K (k) is the convergence power spectrum, 
which depends on the projected 3-dimensional mass power 
spectrum P 3 £>(fc): 



5. Analysis of the systematics 

Now we have to check that the known systematics cannot 
be responsible for the signal. In the following we discuss 
three types of systematics: 

• The intrinsic alignment of galaxies which could exists 
in addition to the lensing effect. We assume such an 
alignment do not exist, but the overlapping isophotes 
of close galaxies produces it. We could in principle re- 
move this effect by choosing a window function small 
compared to the galaxy distance in the pair, such that 
close galaxies do not influence the second moment 
calculation of themselves. However this is difficult to 
achieve in practice. 

• The strongest known systematics is the PSF 
anisotropy caused by telescope tracking errors, the op- 
tical distortion, or any imaginable source of anisotropy 
of the star ellipticity. We have to be sure that the PSF 
correction outlined in Section 2.2 removes any correla- 
tion between galaxy and star ellipticities. 

• The spurious alignment of galaxies along the CCD 
frame lines/columns. We cannot reject this possibil- 
ity since charge transfer along the readout directions 
is done by moving the charges from one pixel to the 
next pixel and so forth, with an transfer efficiency of 
0.99998. This effect could spread the charges of the 
bright objects (very bright and saturated stars pro- 
duce this kind of alignment, but they have been re- 
moved during the masking procedure). Therefore we 
can expect the objects to be elongated along the read- 
out direction. 



P K (k) = \nl f H 4^-P 3D (j^riw) X 

4 J a 2 {w) \Jk{w) J 



Jk(w) is the comoving angular diameter distance out to 
a distance w (wh is the horizon distance), and n(w(z)) 
is the redshift distribution of the sources. The nonlinear 
mass power spectrum Pzo{k) is calculated using a fitting 
formula (Peacock & Dodds 1996). 

We see in Figure 1 that the measured signal is con- 
sistent with the theoretical prediction, both in amplitude 
and in shape. In order to have a better idea of how sig- 
nificant the signal is we can compare for each smoothing 
scale the histogram of the shear variance in the random- 
ized samples and the measured signal. This is is shown in 
Figure 2, for all the smoothing scales shown in Figure 1. 
The signal is significant up to a level of 5.5(7 . Note that 
the measurement points at different scales are correlated, 
and that an estimate of the overall significance of our sig- 
nal would require the computation of the noise correlation 
matrix between the various scales. 

7 with a source redshift distribution which peaks at 0.9 



5.1. Systematics due to overlapping isophotes 

Let us consider the first point in the above list of system- 
atics. In order to study the effect of close galaxy pairs, we 
measured the signal by removing close pairs by varying a 
cut-off applied on the respective distance of close galax- 
ies. Figure 3 shows the signal measured when successively 
closer pairs with d = (no pair rejection), 5, 10 and 20 
pixels have been rejected. The cases d = and d — 5 show 
an excess of power at small scales compared to d = 10 and 
d = 20 (the latter two give the same signal). Therefore we 
assume that for d > 10 we have suppressed the overlap- 
ping isophotc problem, and in the following we keep the 
d = 10 distance cut-off, which gives us a total of ~ 168000 
galaxies for the whole data sets, as already indicated at 
the end of Section 3.2. By removing close pairs of galax- 
ies, we also remove the effect of possible alignment of the 
ellipticities of galaxies in a group due to tidal forces. 

5.2. Systematics due to the anisotropic PSF correction 

We next study the second point concerning the residual 
of the PSF correction. Figure 11 shows that the star ellip- 
ticity correction is efficient in removing PSF anisotropics. 
The raw star ellipticity can be as large as 20% in the 
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Fig. 3. The thick solid line shows the signal as plotted on 
Figure 1. It was obtained with a catalogue of galaxies where 
galaxies closer than 10 pixels were rejected. The three other 
curves show the same signal measured with different rejection 
criteria: the diamond-dotted line is for no rejection at all, the 
triangle-dashed line for galaxies closer than 5 pixels rejected, 
and the square-dashed line for galaxies closer than 20 pixels 
rejected. This figure illustrates that the overlapping isophotes 
of close galaxies tends to overestimate the shear. 



Object size 

Fig. 5. Average galaxy ellipticity (ei) (solid line) and (e2) 
(dashed line) as a function of the object size rn- It is shown 
that the systematic bias of —1% along the ei component is 
fairly galaxy independent. 
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Fig. 4. Average galaxy ellipticity (e a ) versus the average star 
ellipticity {e* a ) for both components a = 1,2. The dashed lines 
are obtained from the fully corrected galaxy ellipticities, as 
given by Eq.(5). The solid lines are obtained from the galaxy el- 
lipticities corrected from the seeing, but without the anisotropy 
correction term P sm p of Eq.(5). Each ellipticity bin contains 
about N = 16000 galaxies, and the error bars are calculated 
assuming Gaussian errors oc N. Except for a constant tiny bias 
along the ei direction, the corrected galaxies are uncorrelated 
with the stellar ellipticity, which demonstrates that the PSF 
correction method works well. 



Fig. 6. Average tangential galaxy ellipticity (74) and radial 
galaxy ellipticity (7,,) versus the distance from the optical cen- 
ter r. As for Figure 4, the dashed lines are obtained from 
the fully corrected galaxy ellipticities, as given by Eq.(5), and 
the solid lines are obtained from galaxy ellipticities corrected 
from seeing but where the anisotropy correction term P sm p 
has been removed from Eq.(5). Each ellipticity bin contains 
about N = 24000 galaxies, and the error bars are calculated 
assuming Gaussian errors oc N. The absence of a significant 
amplitude between the dashed and the solid lines show that 
the optical distortion effect is a negligible contribution to the 
PSF anisotropy. 
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most extreme cases. Figures 12 to 19 show maps of the 
uncorrected and the corrected star ellipticities. The same 
camera used at different times clearly demonstrates that 
the PSF structure can vary a lot in both amplitude and 
orientation, and that it is not dominated by the optical 
distortion (as we can see from the location of the optical 
center, given by the dashed cross). Individual CCD's are 
2K x AK chips, easily identified by the discontinuities in 
the stellar cllipticity fields. 

Next, let us sort the galaxies according to the increas-, 
ing stellar ellipticity, and bin this galaxy catalogue such 
that each bin contains a large number of galaxies. We 
then measure, for each galaxy bin, two different averaged 
galaxy ellipticities (e a ): one is given by Eq.(5) and the 
other by Eq.(5), without the anisotropy correction P sm p. 
The former should be uncorrelated with the star ellipticity 
if the PSF correction is correct (let us call (e a ) this av- 
erage); and the latter should be strongly correlated with' 
the star ellipticity, (let us call (e™ 1 ) this average). Since 
the galaxies are binned according to the stellar elliptic- 
ity, galaxies of a given bin are taken from everywhere in 
the survey, therefore the cosmic shear signal should van- 
ish, and the remaining possible non vanishing value for 
(ei) and (62) should be attributed to a residual of star 
anisotropy. Figure 4 shows (ei) and fa) (dashed lines) and 



anisotropy clearly does not follow the optical distortion 
pattern). 

5.3. Systematics due to the CCD frames 



-1 



) and 



ani\ 



(solid lines) versus respectively (ef 



and (e| tars ). The solid lines exhibit a direct correlation 
between the galaxy and the star ellipticities, showing that 
the PSF anisotropy does indeed induce a strong spuri- 
ous anisotropy in the galaxy shapes of a few percents. 
However, the dashed lines show that the corrected galaxy 
ellipticities are no longer correlated with the star elliptic- 
ity, the average (ei) fluctuates around —1%, while fa) is 
consistent with zero. This figure shows the remarkable ac- 
curacy of the PSF correction method given in KSB. Error 
bars in these plots are calculated assuming Gaussian er- 
rors for the galaxies in a given bin. The significant offset 
of (ei) of 1% might be interpreted as a systematic induced 
by the CCD, as we will see in the next Section, and can be 
easily corrected for. Figure 5 shows that this systematic is 
nearly galaxy independent, and affect all galaxies in the 
same way. This is also in favor of the CCD-induced sys- 
tematic, since we expect that a PSF-induced systematic 
(which is a convolution) would depend on the galaxy size. 

Figure 6 shows the same kind of analysis, but instead 
of sorting the galaxies according to the star ellipticity am- 
plitude, galaxies are now sorted according to the distance 
r from the optical center. The average quantities we mea- 
sure are no longer (ei) and fa) versus (ef ars ) and (e| tars ), 
but the tangential and the radial ellipticity fa) and (e r ) 
versus r. This new average is powerful for extracting any 
systematic associated with the optical distortion. Figure 
6 shows that the systematics caused by the optical distor- 
tion are a negligible part of the anisotropy of the PSF, as 
we should expect from Figures 12 to 19 (where the PSF 
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Fig. 7. Average galaxy ellipticity {e a ) versus the X and Y 
location on the CCDs. As for Figure 4, the dashed lines are 
obtained from the fully corrected galaxy ellipticities, as given 
by Eq.(5), and the solid lines are obtained from the corrected 
galaxy ellipticities where the anisotropy correction term P sm p 
has been removed from Eq.(5). The systematic negative mean 
value of (ei) along lines or columns of the CCD (the two left 
panels) show that the galaxies are preferentially aligned with 
the columns of the CCD in the whole survey. A positive sys- 
tematic value for fa) (the two right panels) is also visible, 
although much less significant. 



Using the same method as in the previous Section, we 
can also investigate the systematics associated with the 
CCD line/columns orientations. Here, instead of sorting 
the galaxies according to the star ellipticity or the distance 
from the optical center, the galaxies are sorted according 
to their X or Y location on each CCD frame. By averaging 
the galaxy ellipticities (ei) and fa) in cither X or Y bins, 
we also suppress the cosmic shear signal and keep only 
the systematics associated with the CCD frame. Figure 7 
shows (ei) and fa) (dashed lines) and (ei) and fa) (solid 
lines) versus (X) and (Y). The plots from the top- left to 
bottom- right correspond respectively to (ei) versus (X), 
fa) versus (A), (ei) versus (Y), and fa) versus (Y). We 
see that (ei) is systematically negative by ~ —1% for both 
X and Y binnings, while fa) does not show any significant 
deviation from zero. This result is fully consistent with the 
dashed lines in Figure 4 which demonstrate that the —1% 
systematic is probably a constant systematic which affects 
all the galaxies in the same way, and which is not related to 
the star anisotropy correction. The origin of this constant 
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shift is still not clear, it might have been produced during 
the readout process, since a negative (ei) corresponds to 
an anisotropy along columns of the CCDs. 

5.4- Test of the systematics residuals 

The correction of the constant shift of —1% along (ei) 
has been applied to the galaxy catalogue from the begin- 
ning. It ensures that there is no more significant residual 
systematic (cither star anisotropy or optical distortion or 
CCD frame), and demonstrates that the average level of 
residual systematics is small and much below the signal. 
However we have to check that the systematics do not 
oscillate strongly around this small value. If it were the 
case, then this small level of systematics could still con- 
tribute significantly to the variance of the shear. This can 
be tested by calculating the variance of the shear in bins 
much smaller than those used in Figure 4 to calculate the 
average level of residual systematics. In order to decide 
how small the bins should be we can use the number of 
galaxies available in the apertures used to measure the sig- 
nal, for a given smoothing scale. For example for 6 = 0.66' 
there is 45 galaxies in average, and for 9 — 3.3' there is 
1100 galaxies. We can therefore translate a bin size into a 
smoothing scale, via the mean number of galaxies in the 
aperture. We found that the variance of the shear mea- 
sured in these smaller bins is still negligible with respect 
to the signal, as shown by Figure 8. The three panels from 
top to bottom show respectively the star anisotropy case, 
the optical distortion case and the CCD frame case. On 
each panel, the thick solid line is the signal with its error 
bars derived from 1000 randomizations. The short dashed 
lines show the ±la of these error bars centered on zero. 
On the top panel the two thin solid lines show (j 2 ) re- 
spectively measured with the galaxies sorted according to 
e stars an j to e stars rpj^ so \id \[ ne [ n the middle panel 

shows {yf } measured from the galaxies, sorted according 
to their distance from the optical center, and the two thin 
solid lines in the bottom panel show (-f 2 ) measured on the 
galaxies sorted according to X and Y. 

In all the cases, the thin solid lines are consistent with 
the ±lcr fluctuation, without showing a significant ten- 
dency for a positive (j 2 ). We conclude that the residual 
systematics are unable to explain the measured (j 2 ) in 
our survey, and therefore our signal is likely to be of cos- 
mologicai origin. 

A direct test of its cosmologicai origin is to measure 
the correlation functions (e t (0)e t (6>)}, (e r (O)e r (0)) and 
(e r (0)et(0)), where et and e r are the tangential and ra- 
dial component of the shear respectively: 

e t = -ei cos(20 fe ) - e 2 sin(26> fc ) 

e r = -e-2 cos(26> fe ) + ei sin(2# fe ), (16) 

where 9k is the position angle of a galaxy. If the signal is 
due to gravitational shear, we can show (Kaiser 1992) that 



(et(0)e t (8)) should be positive, (e r (O)e r (0)) should show 
a sign inversion at intermediate scales, and (e r (O)et(0)) 
should be zero. This is a consequence of the scalar ori- 
gin of the gravitational lensing effect and of the fact that 
galaxy ellipticity components are uncorrelated. Although 
we do not yet have enough data to perform an accurate 
measurement of these correlation functions, it is inter- 
esting to check their general behavior. Figure 9 shows 
that in our data set, although the measurement is very 
noisy, both (e t (O)e t (0)) and (e r (O)e r (0)) are positive val- 
ued, while (e r (O)e t (0)} is consistent with zero. This mea- 
surement demonstrates that the component of the galaxy 
ellipticities e a of well separated galaxies are uncorrelated, 
and it is in some sense a strong indication that our signal 
at small scales is of cosmologicai origin. 

The last thing we have checked is the stability of the 
results with respect to the field selection. We verified that 
removing one of the fields consecutively for all the fields 
(see Section 2 for the list of the fields) docs not change 
the amplitude and the shape of the signal, even for the 
Abell 1942 field. The cluster has no impact and does not 
bias the analysis because it was significantly offset from 
the optical axis. This ensures that the signal is not pro- 
duced by one field only, and that they are all equivalents in 
terms of image quality, PSF correction accuracy and sig- 
nal amplitude, even using V and I colors. It also validates 
the different pre-reduction methods used for the different 
fields. 

6. Cosmologicai constraints 

Figure 1 provides a first comparison of our signal with 
some cosmologicai models. In order to rule out models we 
need to estimate first the sample variance in the variance 
of the shear. Although it has not been yet exactly de- 
rived analytically (because calculations in the non-linear 
regime are difficult), ray-tracing simulations can give an 
accurate estimate of it. We used the ray-tracing simula- 
tions of Jain et al. 2000 for this purpose. 



Table 2. List of the ray tracing simulations we used (see 
Jain et al. 2000 for details). The redshift of the sources is 1. 



Simulation ^ 


r 


fio 


A 


erg 


(1) OCDM 


0.21 


0.3 





0.85 


(2) rCDM 


0.21 


1 





0.6 


(3) rCDM 


0.21 


1 





1 



Table 2 shows the two simulations we used. The rCDM 
model with a% = 1 is not an independent simulation, but 
was constructed from the rCDM model with as = 0.6 sim- 
ply by dividing k by 0.6. This should empirically mimic 
a model with both f2 and cr 8 equal to one. The redshift 
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Fig. 8. Possible contribution of the systematics studied in Sec- 
tion 3.2 to the signal. On each of the plots, the thick solid line 
shows the signal as displayed on Figure 1 , and the dashed lines 
show the ±1(7 fluctuation obtained from 1000 random realiza- 
tions. From top to bottom: (a) The two thin solid lines are 
(7 2 ) 1 / 2 measured on the galaxies sorted according to the star 
ellipticity strength (see Figure 4) . For the different smoothing 
scales, the mean number and the variance of the number of 
galaxies in the chosen bins fit the one observed in the signal 
(thick solid) curve, (b) the thin solid line is (7?) 1 ^ 2 measured 
on the galaxies sorted according to the distance from the op- 
tical center, and on (c) the two thin solid lines correspond to 
the galaxies sorted according to their X or Y location on the 
CCDs. 
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Fig. 9. From top to bottom, measurement of the correlation 
functions (e t (O)e t (0)), (e r (O)e r (0)) and (e r (O)e t (0)). The error 
bars are computed from 50 random realizations of our data set 
where the orientations of the galaxies were randomized. 



of the sources is equal to 1, which is not appropriate for 
our data. However, for the depth of the survey, we be- 
lieve that it represents fairly well the mean rcdshift of 
the galaxies, which is the dominant factor in determining 
the second moment. Figure 10 shows the amplitude and 
the scale dependence of the variance of the shear for the 
three cosmological models, compared to our signal. It is 
remarkable that models (1) and (3) can be marginally re- 
jected (We did not plot the error bars due to the intrinsic 
ellipticity for clarity: they can be obtained from Figure 8). 
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Our measurements are in agreement with the clus- 
ter normalized model (2). Also plotted is the theoretical 
prediction of a ACDM model, with Cl = 0.3, A = 0.7, 
r = 0.5 and a redshift of the sources z s = 1. It shows 
that the low-Jl model is also in good agreement with the 
data, which means that weak gravitational lensing pro- 
vides cosmological constraints similar to the cluster abun- 
dance results (Eke ct al. 1996,Blanchard et al. 1999): the 
second moment of the shear measures a combination of 
as and fio (see equation 8). A measure of the third mo- 
ment of the convergence would break the fl-as degeneracy, 
but this requires more data (see Bernardeau et al. 1997, 
Van Waerbeke et al. 1999, Jain et al. 2000). It should also 
be noted that for the simulations, we have considered 
cold dark matter models with shape parameter T = 0.21; 
higher values of T increase the theoretical predictions on 
scales of interest, e.g. the fio = 1, <7s = 1 model would be 
ruled out even more strongly. We conclude that our anal- 
ysis is consistent with the current favored cosmological 
models, although we cannot yet reject other models with 
high significance. Since we have only analyzed 2 square 
degrees of the survey, with forthcoming larger surveys we 
should be able to set strong constraints on the cosmolog- 
ical models as discussed below. 

Due to the imprecise knowledge of the redshift distri- 
bution in our data, the interpretation might still be sub- 
ject to modifications. The final state of our survey in 4 
colors will however permit the measurement of this distri- 
bution by estimating photometric redshifts for the source 
galaxies. 

7. Conclusion 

We have demonstrated the existence of a significant cor- 
relation between galaxy ellipticities from 0.5 to 3.5 arc- 
minutes scales. The signal has the amplitude and the an- 
gular dependence expected from theoretical predictions 
of weak lensing produced by large-scale structures in the 
universe. We have tested the possible contribution of sys- 
tematic errors to the measured signal; in particular we 
discussed three potential sources of spurious alignment of 
galaxies: overlapping isophotes of very close galaxies, star 
anisotropy and CCD line/column alignment. The first of 
these systematics is easy to deal with, simply by remov- 
ing close pairs, although we may have decreased the signal 
slightly by removing them. The star anisotropy seems to 
be very well controlled, in part due to the fact that the 
bias adds quadratically with the signal. Moreover, in the 
absolute sense, the bias does not exceed a fraction of 1 
percent, which is adequate to accurately measure a vari- 
ance of the shear of few percent. The only important bias 
we found seems to be associated with the CCD columns, 
and it is constant over the survey, it is therefore easy to 
correct for. The origin of this CCD bias is still unclear. 

As an objective test of the reality of the gravitational 
shear signal, we measured the ellipticity correlation func- 
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Fig. 10. Comparison of our signal (thick line) with three cos- 
mological models. The error bars are the cosmic variance mea- 
sured on five independent realizations at the smoothing scale 
indicated by the x-axis. For clarity, the shot noise error bars 
of the signal are not plotted, their amplitude can be read on 
Figure 8. From bottom to top, the dashed lines correspond to: 
model (1), model (2) and model (3) as given in Table 2. The 
shot-noise error bars of the signal are in fact comparable in am- 
plitude to the cosmic variance error bars of model (2). We show 
also a cluster-normalized A model (dotted line) with £1 = 0.3, 
A = 0.7, and a CDM power spectrum with T = 0.5. This model 
was not obtained from a simulation, but computed using the 
non-linear power spectrum using the Peacock & Dodds 1996 
formula. 



tions (e t (O)e t (0)>, (e r (O)e r (0)) and (e r (O)e t (0)). While the 
measurement is noisy, the general behavior is fully consis- 
tent with the lensing origin of the signal. The tests for sys- 
tematic errors and the three ellipticity correlation function 
measurements described above have led us to conclude 
with confidence that we have measured a cosmic shear 
signal. 

With larger survey area, we expect to be able to mea- 
sure other lensing statistics, like the aperture mass statis- 
tic (M ap ; see Schneider et al. 1998). The M ap statistic is 
still very noisy for our survey size because its signal-to- 
noise is lower than the top-hat smoothing statistic, due to 
higher sample variance (We verified this statement using 
the ray tracing simulation data of Jain et al. 2000). Our 
survey will increase in size in the near future (quickly up to 
7 square degrees) , leading to a factor of 2 improvement in 
the signal-to-noise of the results presented here. According 
to our estimates, this will be enough to measure M ap at 
the arc-minute scale with a signal-to-noise of ~ 3. The de- 
tection of the skewness of the convergence should also be 
possible with the increased survey area. This will be im- 
portant in breaking the degeneracy between the amplitude 
of the power spectrum and (Bernardeau et al. 1997, 
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Van Waerbeke et al. 1999, Jain et al. 2000). These mea- 
sures should also provide nearly independent confirma- 
tions of the weak gravitational lensing effect as well as 
additional constraints on cosmology. Thus by combin- 
ing different measures of lensing by large-scale struc- 
ture (top-hat smoothing statistics, M ap statistics, cor- 
relation function analysis, power spectrum measure- 
ments), higher order moments, and peak statistics 
(Jain & Van Waerbeke 2000), from forthcoming survey 
data, we hope to make significant progress in measuring 
dark matter clustering and cosmological parameters with 
weak lensing. 

We also hope to do a detailed analysis with a more 
sophisticated PSF correction algorithm. For instance, the 
mass reconstruction is linear with the amplitude of the 
residual bias, and a fraction of percent bias is still enough 
to prevent a definitive detection of filaments or to map 
the details of large scale structures. Since we show else- 
where (Erben et al. 2000) that such a bias is unavoid- 
able with the present day correction techniques and image 
quality, there is still room to improve the analysis prior 
to obtaining accurate large-scale mass maps. Recent ef- 
forts to improve the PSF correction are very encouraging 
(Kaiser 1999). We plan to explore such approaches once 
we get an essentially homogeneous data set on a larger 
field. 
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Fig. 11. Star ellipticities of all the survey before (top panel) and after (bottom panel) the correction. After correction, the star 
ellipticity is randomly distributed around zero, as expected. 




Fig. 12. Uncorrected (left) and corrected (right) star ellipticities for FIELD F14P1. The dashed cross shows the location of the 
optical center. Frames are graduated in pixels. The reference stick at the top-left of the frame shows the amplitude of a 10% 
distortion. This length of reference applies also for Figures 13 to 19. 
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Fig. 15. Same as Figure 12 for FIELD 03hrie. 
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Fig. 16. Same as Figure 12 for FIELD SA57. 
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Fig. 17. Same as Figure 12 for FIELD al942. 
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Fig. 18. Same as Figure 12 for FIELD F02P1. 
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Fig. 19. Same as Figure 12 for FIELD F02P4. 
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